Majorana fermion chain at the Quantum Spin Hall edge 
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We study a realization of a Id chain of Majorana bound states at the interfaces between ahernating 
ferromagnetic and superconducting regions at a quantum spin HaU insulator edge. In the hmit of weU 
separated Majoranas, the system can be mapped to the transverse field Ising model. The disordered 
critical point can be reached by tuning the relative magnitude or phases of the ferromagnetic and 
superconducting order parameters. We compute the voltage dependence of the tunneling current 
from a metallic tip into the Majorana chain as a direct probe of the random critical state. 
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Some electronic phases support emergent quasiparti- 
cle excitations that are in a sense fractions of the orig- 
inal electrons, as in the fractional quantum Hall effect. 
One of the most basic fractional excitations is the Ma- 
jorana fermion familiar from particle physics: the Majo- 
rana fermion is its own antiparticle and represents "half" 
of an ordinary (Dirac) fermion. Several condensed mat- 
ter systems arc theoretically believed to realize Majorana 
fermions in their excitation spectrum, and their direct 
observation is a major goal of current research. In this 
work we compute the tunneling conductance and other 
experimental signatures of one of the simplest systems of 
many Majorana fermions. 

Proposals to observe Majorana fermions [1-3] have 
been based on the proximity effect between an ordinary s- 
wave superconductor and the recently discovered "topo- 
logical insulator" phases in two-dimensional [4-6] and 
three-dimensional [7-10] materials. These phases gen- 
erally support gapless edge or surface states. In the two- 
dimensional (2d) case, the gapless edge can be viewed 
as countcrpropagating up-spin and down-spin electrons, 
hence the terminology "quantum spin HaU" (QSH) for 
this phase. In the three-dimensional "strong topologi- 
cal insulator" [9], the surface state is a two-dimensional 
metal with a Fermi surface that encloses an odd num- 
ber of Dirac points. The proximity effect between an 
ordinary superconductor and this surface state leads to a 
superconducting state [2] that is time-reversal symmetric 
but topologically similar to the p+ip superconductor [11] , 
with Majorana fermions trapped in vortices. 

Perhaps the most appealing from a technical perspec- 
tive are Majorana realizations in Id systems. In partic- 
ular, a QSH edge supports a localized Majorana fermion 
excitation [1, 12]. When the edge state is gapped by 
coupling to a superconductor (SC) in one region and a 
ferromagnct (FM) in another, a single zero-energy Ma- 
jorana fermion appears at the SC/FM boundary if it is 
sufficiently narrow (Fig. 1(a)). In this paper we study 
how unique features of a chain of Majorana fermions, 
created by alternately FM and SC regions along a QSH 
edge, can be observed using a metallic tunneling tip. 
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FIG. 1: (a) A Majorana chain realization: An array of al- 
ternating FM and SC regions at the QSH edge. Following a 
real-space RG analysis, the chain decouples into pairs, shown 
by the upper lines, is the energy of the pair decimated with 
the i-th highest energy, (b) The RG effectively couples the 
STM tip to every Majorana pair, with tunneling suppressed 
by the coupling coefficient a. Majorana pairs decimated at 
the edge are coupled more strongly to the lead than pairs 
decimated in the bulk of the chain. 



An ideal, uniform, system realizes the interesting free- 
Ma jorana fermion quantum critical point. Realistically, 
however, FM and SC domains will inevitably vary ran- 
domly in size and proximity amplitude, and therefore 
the proposed system would realize an even more inter- 
esting model: the random Majorana-fermion chain in its 
random-singlet phase [13]. Previously, likely experimen- 
tal Majorana signatures were discussed for one or two 
Majorana pairs [3, 14-16], and the phase diagram of a 
2d Majorana fermion model without disorder was ob- 
tained [17]. Experimentally, it may be even simpler to 
pattern an irregular phase-separated mixture of SC and 
FM materials to contact the QSH edge, rather than a 
precise controlled-lithography FM-SC-FM configuration. 

Furthermore, the Majorana fermion chain may give a 
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first convincing measurement of tlie random-singlet phase 
and its associated Griffiths scaling. Contrary to spin 
chains where neutron scattering or other spin-sensitive 
probes can only give limited information, the Majorana 
chain setup allows a direct electronic tunneling measure- 
ment that reveals the nature of the random singlet phase. 
In what follows we will first derive the Hamiltonian of the 
random Majorana chain contacted by a tunneling tip, 
and then calculate its main tunneling I-V features. 

The Majorana chain Hamiltonian arises as the low en- 
ergy (in-gap) sector of the QSH edge Hamiltonian in 
proximity to SC and FM domains alternately, with in- 
terfaces arising at xi = 0, and Xi > thereafter. The 
QSH edge Hamiltonian has the form [2] H = i^'^'H^', 



where ^! 



and 



n = -ivd^a'^T'' + A(x) ■ f + M{x) ■ a 



(1) 



The first term describes the counterpropagating spins of 
the QSH edge, while the rest describes the proximity 
to FM regions {i = 2j) and SC regions {i = 2j + 1) 
with M{x) = {Mq cos (j)mi, Mq sin (j)mi,0) and A(a;) = 
{Aq cos (j) si, Aq sin (p si, 0)- a^'^'^ are Pauli matrices that 
act in the space of right and left movers V't and ipi, while 
^x,y,z g^j.g pg^^jj matrices that act on the ip and V'^ blocks. 
For simplicity, we set all magnetic moments Mq and pair- 
ing strengths Aq uniform, but let the x-y moment direc- 
tions (j)M2j and SC order parameter phase 0s2j+i vary. 
Magnetic moment along z does not open a gap and is ir- 
relevant for what follows. Each interface at Xi will host a 
Majorana state, whose creation operator we will denote 
7i. Unlike standard fermion creation operators, a Majo- 
rana fermion obeys j'^ = j and 7^ = 1. Two Majorana 
fermions can combine to give a Dirac fermionic state via 



'H12 = -«ii27i72 ~ ti2(2c^c - 1) 



(2) 



In a Majorana chain such an interaction weakly couples 
nearest neighbor states. 

It is instructive to derive this interaction for just two 
Majoranas located at a;i = and X2 = L due to a FM- 
SC-FM sequence, described by (/)j\/o, 0si, 0m2 respec- 
tively. To obtain the Majorana state at xi = we ignore 
the X2 interface and find a zero energy Majorana solution 
of the Hamiltonian (1) localized about a; = 0: 
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exchange A and M in Eq. (1) switches the second and 
third components ^(2) ^(3) in the resulting wave func- 
tion vector. Thus the X2 = L Majorana is: 
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with = v/Mq for x > L, and £- ~ v/Aq for x < L. 
Once again we define 72 = / dx^2{x)'^{x). 

The Majorana wave function ^1 is a zero-energy eigen- 
state of (1) ignoring all interfaces other than xi. The 
interface at X2 clearly gives rise to hybridization of 71 
with 72 . This we obtain by writing H = Ho + V with Hq 
describing just the FM/SC interface at xi and: 



V^eix- L)(Ai • f,y - M2 ■ T.,y). 



(5) 



Q{x) is the Heaviside function. The perturbation piece 
V produces the hybridization 



^16) = / dx^lvu^). 



(6) 



Writing this hybridization so that it captures SC/FM/SC 
as well as FM/SC/FM Majorana pairs, we obtain 



tij+i = 2 cos 



A0, 



Ao + Mo 



(7) 



the transformation c = ^(71 + 172), c'^ = ^(71 — 172); a Int 
pair-wise Majorana interaction is simply written as: 



with A(j)2j+i = (j)M2j+2 - 4>M2] and £ = v/Aq de- 
scribing a FM/SC/FM sequence, and A(/)2j = <j>s2j-i — 
0S2J+1 and e = v/Mo for a SC/FM/SC sequence. 
Since the interaction between Majorana fermions de- 
cays exponentially with distance, we keep only the 
nearest neighbor interaction, and the Hamiltonian of 
the chain becomes ^ . ti.i_|_i7i7i-|_i. The strong-disorder 
random-singlet fixed point is attained when lnt2j,2j+i = 



2j-|-l,2j+2 



[18] 



Mapping this Hamiltonian to the quantum Ising spin 
chain [13] immediately yields critical thermodynamic [18] 
and entanglement [19] properties of the Majorana chain, 
which seem difficult to measure in practice. Instead we 
focus on the local density of states, which is probed di- 
rectly using a single STM lead biased at voltage V rel- 
ative to the chain. The combined Hamiltonian for the 
Majorana chain and a probing tip coupled to one end is: 



H = -i 



1=1 
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-lf(3) 



with the decay lengths for x > and x < being 
= v/Ao and £- = v/Mq respectively. The Majo- 
rana creation operator is thus 71 = / dx^i{x)'^{x). Du- 
ality yields the second Majorana state at X2 = L: an 
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where ti i^i is given in Eq. 7, and tJj is the annihilation 
operator of tip electrons [14] . We will use real-space RG 
to reduce the semi-infinite chain and STM to aproblcm 
of a tip interacting with independent Majorana pairs: 

00 

n=l 

(9) 
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The decoupled Majorana pairs correspond to the so- 
called random-singlets formation that was investigated 
in random Heisenberg chains [20] and in random hop- 
ping fermions [21]. The I-V characteristics of this Hamil- 
tonian are clear. Whenever the voltage of the tip is at 
resonance with ±e„ a step appears in the I V curve, 
with an amplitude proportional to the tunneling rate, a^: 



dl 
dV 



aliSieV + €„)-SieV -e„)). 



(10) 



The two opposite sign delta-functions are due to the spe- 
cial coupling term oc "0 -|- V-"^ ; if at a positive resonance, 
y = ]e„], i/; is responsible for a current flowing into the 
chain, and when V = — jenj, ip^ will produce the current, 
which flows out o/thc chain to the tip. 

The measured LDOS, Eq. (10), will provide not only 
the chain's density of states, but also their overlaps a„ 
with the edge Majorana. Therefore it reveals both the en- 
ergy and spatial structure of the random Majorana chain. 
We can find both the energy levels e„ and the edge over- 
lap, a„ using the real space RG method applied 
series of unitary transformations on Eq. 8 [22]. Its ad- 
vantage is that it does not eliminate pieces of the Hilbert 
space, but rather decomposes the Hamiltonian Eq. (8) 
into decoupled Majorana pairs. 

The unitary real-space RG procedure finds a sequence 
of unitary transformations which isolate the strongest 
bonds, at each stage, in Eq. (8) from the rest of the 
chain. For instance, if bond n is the strongest, we find 
that using S^") = ^f^ijn-iln+i - yields: 



* 57 7n-l7n+2 



(11) 

to order O ^^-^ . Thus, bond n is decoupled, and bridg- 
ing over it, a new coupling between sites n — I and n + 2 
appears. Iteratively applying the transformations Si, 
which decouple the sites Uia and n^fc, also yields a wave 
function for the interacting chain, which is just a trans- 
formed product wave function of the decoupled pairs: 



(12) 



To find the LDOS, wc need to compute ground state 
matrix elements of the contact term in Eq. (8). Our 
method allows us to write such matrix elements in 
terms of the decoupled Majorana-pairs wave functions, 
|'I'„.^„.j,), as long as we use the transformed contact term. 
Applied to the contact term, the unitaries turn 71 into a 
linear combination of the other Majorana's in the chain: 



-iSi 



anlleadln (13) 
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The bipartite geometry of the chain guarantees that only 
odd sites can couple to the lead. At the end of the RG 



process, the STM lead effectively couples to a Majorana 
site in many decoupled Majorana pairs, and the total cur- 
rent is the sum of the current into each pair (Eq. 10). The 
dominant contributions, however, arise when the active 
edge of the chain (at any stage of the RG) are decimated. 
E.g., consider the leftmost Majorana, 71. Its decimation 
due to a strong bond with 72, leads to 

t2 



e'^-fie ^ 71 
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(14) 



Subsequent decimations of the new edge 73 couples the 
tip to the next odd-numbered Majorana site which is the 
effective left edge at that energy scale. Bulk decimations 
not containing the edge produce only subdominant corre- 
lations contributions in Eq. (13). The tunneling dl/dV 
will have strong spikes corresponding to the edge decima- 
tions, between which, bulk decimations yield additional 
spikes of various strengths (but always weaker than both 
edge spikes; see below). 

Next, we find the probability distribution for the cou- 
pling coefficients a„, and the energies of the Majorana- 
pairs they connect to the edge. These are obtained 
from the universal coupling distributions of the effec- 
tive bond strength at low energies. Define — max{ti} 
as the maximum coupling in the chain at a given point 
in the RG process. Also, define the logarithmic bond 
variables: Q = hi ft /ti, and logarithmic fiow parameter 
r = \nili/il, with ilj the largest initial coupling. Af- 
ter many bond decimations, the distribution of bonds 
settles into the universal random-singlct fixed point dis- 
tribution: [18, 20] 



PiC) 



(15) 



We define u„(t; F) as the probability that when the RG 
fiow parameter reaches F, the active edge Majorana is 
coupled with amplitude e~'^ to the lead, and is the result 
of n previous decimations involving the edge Majorana, 
as outlined in Eq. (14). The evolution equation for u„ is 
given by = 

-P0U„(t)+ J dT'dC2PoUn-l{r')P{C2)S{T-T'-C2), (16) 

where Pq = P{0), and the probability of a bond to be 
decimated as F changes by dT is PodT. The first term in 
Eq. 16 marks the decrease in this probability as the edge 
is decimated. The second is a source term for the n'th 
edge, passing the correlation information from the edge 
after n — 1 decimations. The delta function expresses the 
reduction in correlation of the new edge compared to its 
predecessor, by a factor e~^^ = t2/^ [as in Eq. (14)]. 
il = ti since the first bond is decimated, and a factor of 
2 was neglected since to logarithmic accuracy we expect 
C2 > ln(2) (see, e.g., [23]). The solution to Eq. 16 is 



Fo 1 „ /A + l/Fo 

Uji = r In ; — 

F n! V ^ + l/r 



(17) 
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FIG. 2: Marginal distribution of edge energies for n — 1-4. 
The n'th peak has a maximum at ln(ln(l/ii')) — n/2. dl/dV 
thus has peaks at regular intervals in ln(ln(l/l/)). 
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FIG. 3: Typical IV curve obtained by integrating dl /dV from 
high voltage down to V. A sharp rise in current occurs at 
values of V corresponding to typical energies of edge pairs. 



where A is the Laplace transform variable of r. 

From Un it is easy to state the probability density for 
the n'th edge to be decimated at F and with correlation 
e""^ to the edge. It is simply 



1 



d4'"^Se)fdY = PoUn = -U. 
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(18) 



As the voltage is scanned from the band-edge towards 
the band center, the edge decimation will produce dl/dV 
peaks at values of Vn with a distribution in F„ = 
InQi/Vn given by dA^"^^"^ /dT at A = 0. Finding the 
maximum of this distribution yields the most probable 
voltage of the n'th edge decimation: 



(19) 



The statistical properties of the amplitude of the n'th 
dl/dV peak, which is w„ = tge"^"^, are also encoded in 
dJn'^^'^) / dr{T, F). For instance, the average height of the 
n'th peak, (e~^'^" " is given by integrating over F in 
Eq. (18) to obtain the transform of the marginal dis- 
tribution for r, and then setting A = 2. This yields, 
roughly, (w„) ^. This average, however, is dom- 
inated by samples with anomalously high correlations. 
More useful is the typical weight of the n'th edge peak, 
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peak is Vri 
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^ given that the energy of the n'th 
This yields the behavior: 
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It can be shown [24] that summing up all currents due to 
bulk decimations occuring between the n'th and n + I'th 
edge peaks gives Wn^^'^'^ w^f'^^'^-'e"^^"'^"''', which allows 
us to safely neglect them in the single-lead setup. 

Fig. 3 shows the predicted current plotted vs. 
ln(ln(l/y)). Each time the voltage crosses an edge peak 
from above, there is a sharp rise in the current as the 
probe couples to a decimated Majorana pair. For higher 
n, i.e., edge decimations at lower energies, the typical 
coupling and current jumps decrease sharply. 

In this manuscript we calculated the distributions of 
the STM tunneling characteristic of an experimentally 



feasible realization of an interacting Majorana chain. 
The observation of the unusual tunneling STM current 
we find would be strong evidence for the critical Majo- 
rana chain. The detailed transport information was ob- 
tained using a unitary transformations version of the Ma- 
Dasgupta RSRG procedure, which allows going beyond 
previous studies of transport in random-singlet type crit- 
ical points, where only averages could be obtained [25]. 
This technique could also be applied to the edge tunnel- 
ing conductance of the random critical point of an ordi- 
nary fermionic bipartite chain. In future work we will 
explore this spectrum in more detail, as well as consider 
2-tip transport through a Majorana chain. 
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